Fluctuation Theorem for Arbitrary Open Quantum Systems 
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Based on the observation that the thermodynamic equilibrium free energy of an open quantum 
system in contact with a thermal environment can be understood as the difference between the free 
energy of the total system and that of the bare environment, the validity of the Crooks theorem 
and of the Jarzynski equality is extended to open quantum systems. No restrictions on the nature 
of the environment or on the strength of the coupling between system and environment need to be 
imposed. This free energy entering the Crooks theorem and the Jarzynski equality is closely related 
to the Hamiltonian of mean force that generalizes the classical statistical mechanical concept of the 
potential of mean force. 
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Since its formulation in 1997 the classical nonequilib- 
rium work relation by Jarzynski [l| (now commonly re- 
ferred to as Jarzynski Equality) 
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kept raising questions and concerns on its range of valid- 
ity and applicability. Here w denotes the work performed 
on a system when some parameters of this system are 
changed according to a prescribed protocol. Initially, the 
system is supposed to be prepared in a thermal equilib- 
rium state at the inverse temperature (3. The brackets 
(•) denote a non-equilibrium average over many repeti- 
tions of this process, running under the same protocol. 
According to the Jarzynski equality the average of the 
exponentiated negative work is independent of the de- 
tails of the protocol and solely determined by the thermal 
equilibrium free energy difference IS.F between the initial 
equilibrium state and a hypothetical equilibrium state at 
the initial temperature and those parameter values that 
are reached at the end of the protocol. In the mentioned 
paper [l| the validity of this equality was demonstrated 
within a classical statistical approach for isolated sys- 
tems which initially are in the required equilibrium state 
at inverse temperature f3 and also for classical systems 
that stay in weak contact with a thermal bath during the 
protocol. 

Considerable efforts have been devoted to the devel- 
opment of the quantum version of Eq. ([1]), and more 
generally of the Crooks fluctuation theorem [3 that un- 
derlies it, i.e.. 



Pto,tf{-w) 



(2) 



where ptf_ta{w) denotes the probability density function 
(pdf) of work performed by the parameter changes ac- 
cording to a protocol running between the initial time 
to and final time tf. The pdf of work for the re- 
versed protocol is denoted by ptQ,tf{w). All these at- 
tempts refer to quantum isolated or weakly coupled sys- 



tem with Hamiltonian or Markovian dynamics, respec- 
tively d, i, [1, H, 0, i]. The proof for the validity of Eqs. 
([1]) ; (l2|) in the quantum case with weak coupling allowing 
for an otherwise general non Markovian dynamics of the 
open quantum dynamics and arbitrary force protocols 
was provided only recently in Ref. 0]. 

The applicability of Eqs. ([1]), ([2]) to the case of weak 
coupling is consistent with the construction of quantum 
and classical statistical mechanics which relies on that 
assumption. In striking contrast, extending the meth- 
ods of statistical mechanics to cases that involve a non- 
negligible system-environment interaction presents a ma- 
jor challenge [lol . 11 1. Addressing this question is by now 



becoming more and more pressing, as the advancement of 
technology poses us in the position to investigate exper- 
imentally the thermodynamic behavior of nanosystems 
operating in the quantum regime, whose reduced sizes 
make the system-environment coupling an important is- 
sue. 

A satisfactory generalization of the applicability of 
Eq. ^ for the classical strong coupling regime was 
put forward by Chris Jarzynski himself [12l | ; it should be 
stressed, however, that the objective of the corresponding 
quantum treatment has not been achieved yet. The key 
tool used in Ref. [l3| to overcome the difficulties posed 
by the presence of strong coupling is the Hamiltonian of 
mean force H*{Ts,t) where Tg denotes a point in the 
phase space of the subsystem of interest. This Hamil- 
tonian of mean force is defined as the effective Hamil- 
tonian that describes the Boltzmann-Gibbs equilibrium 
of the marginal probability density of the subsystem of 
interest; it reads 



H*irs;t)^HsiTs,t) 
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1 J dTBeM-l3iHB{TB) + Hsb(Ts,Tb))) 
— in- 
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where P b denotes a point in the phase space of the bath. 
The total Hamiltonian of system plus environment is 
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given by 



quantum Tasaki-Crooks Fluctuation Theorem (sl. [l^: 



H{Ts,TB,t) ^ HsiTs,t) + HBiTB) + HsB{Ts,rB) (4) 

which is composed of the Hamiltonian of the isolated 
system of interest, Hs{Ts,t) (time-dependent), the bath 
Hamiltonian Hb(Xb) and the interaction Hamiltonian 
Hssi^Si^B)- This Hamiltonian of mean force general- 
izes the concept of the potential of mean force, see Eq. 
(|14[1 below, that is commonly employed, for example, in 
reaction rate theory [l3l | and in the study of implicit sol- 
vent models in terms of biomolecular simulations 
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In the context of quantum rate theory, potentials of 
mean force have been determined for the reaction coordi- 
nate from path integral expressions of the partition func- 
tion of a composed system with the reaction coordinate 
confined to a "centroid" 3, 3]. A direct apphcation of 
this very approach to obtain quantum fluctuation theo- 
rems for open systems though is not obvious. Jarzynski, 
in fact, did emphasize that his treatment is restricted to 
the classical case [l2j. Addressing this problem thus re- 
quires a careful analysis of what one should consider as 
the system partition function from which an equilibrium 
free energy of the system can be inferred. One could 
naively take for this partition function the bare parti- 
tion function of the isolated system of interest: This pro- 
cedure would, however, neglect the prominent fact that 
the interaction with the bath alters the system proper- 
ties. Such a choice, therefore, is generally physically not 
correct. Instead of the partition function of the isolated 
system one has to choose a properly defined partition 
function, which embraces the influence of bath on the 
open quantum system. 

As we will show in this Letter, the introduction of a 
proper partition function for an open system allows one 
to prove that both the Tasaki-Crooks fluctuation the- 
orem and the Jarzynski equality in fact hold true for 
general open quantum systems, independent of coupling 
strength. Moreover, from this partition function a quan- 
tum Hamiltonian of mean force can be inferred which 
takes over the role of the Hamiltonian of mean force in 
classical statistical mechanics. 

The Argument.— Consider the Hamiltonian operator 
H{t) of a quantum system composed of the interacting 
system and the bath which we write as 



H{t)^Hs{t)+HsB+HB 



(5) 



where the system Hamiltonian Hs{t) is time-dependent 
in a way that results from a pre-specified protocol of 
system parameter changes. The interaction Hamiltonian 
HsB and the bath Hamiltonian Hb are supposed to be 
independent of time. The change of system's parameters 
can be interpreted as a time dependent external forcing 
that is able to perform work on the system. 

The total system is isolated; it therefore obeys the 



Pto,tf{~w) Y{to) 
with Y{t) being the total partition function, i.e.. 



Y{t) = Tre^0(iisit)+HsB+HB 



(6) 



(7) 



where Tr denotes the trace over the total system Hilbert 
space, and the symbols ptf^toiw) and pta,t]{w) denote 
the probability densities of doing the work w when the 
protocol is run in the forward and backward directions, 
respectively. It is important to note that, due to the 
fact that the forces solely act on the system, the work 
performed on the open system coincides with that done 
on the total system. 

In order to properly define the partition function of the 
open quantum system S staying in thermal equilibrium 
with a bath, we appeal to thermodynamic reasoning. As 



pointed out in Ref. [17| , the free energy of the open sys- 



tem of interest is the difference between the total system 
free energy and the bare bath free energy: 



Fs{t) = F{t) - Fb . 



(8) 



Here t merely specifies the values of the external pa- 
rameters as they occur in the course of the protocol 
at the time t. The function ^^(t) satisfies all re- 
quired properties of a free energy [lOl . 
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120|. Us- 
ing the statistical mechanical relation (3F = —\nZ be- 
tween equilibrium free energy and the partition func- 
tion, one finds from Eq. ([5]) for the partition func- 
tion of an open quantum system the result that is well 
known to those working on strong quantum dissipation 
HI, [13, El, El, 0, [m, M, [23, m , namely: 



Zs{t) = 



rYj.g-PiHsit)+HsB+HB) Y{t) 



Trse 



Zb 



(9) 



where Tr^ is the trace over the bath Hilbert space, 
and Zb is the bare bath partition function, i.e., Zb = 
Tr^e"^^^, which is independent of time. 

From Eqs. ([7]) and (O, and the fact that Zb does 
not depend on time t, the salient relation Y{tf)/Y{to) = 
Zs{tf)l Zs{io) follows. This quantum result assumes a 
form reading just like in in the classical case [l^l. There- 
fore, Eq. ([6]) becomes 



Ptf,to{+w) _ Zsitf) 



n(w-AFs) 



Pto,tf{- 



(10) 



which states that the ratio of probabilities of work in 
the backward and forward protocols is dictated by the 
equilibrium free energy difference of the open quantum 
system, i.e., AFs = Fs{tf) - Fsito), with Fs{t) given by 
Eq. ®. 
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By multiplying Eq. (fTO)) by pto,tf{—w)e~^^ and in- 
tegrating over w in the usual way one obtains with AFs 
the very form ([1]) of the Jarzynski equality for open quan- 
tum systems, being valid independently of the coupling 
strength and the details of the bath. This is in complete 
analogy with Jarzynski's classical result, which, there- 
fore, carries over to the quantum case. Hence, if a clas- 
sical force acts on an open system, the average expo- 
nentiated work e"*^™ equals the exponentiated system 
equilibrium free energy difference, both in classical and 
quantum regimes. 

Remarks.— The system partition function Zs{t), de- 
fined in Eq. is actually the partition function asso- 
ciated to the quantum Hamiltonian of mean force H* (t) , 
defined in analogy to the classical Hamiltonian of mean 
force as: 

H*(t) -.^ In . (11) 

In fact, the partition function Zgit) can be recast as: 



of mean force -, mentioned before: 



Zs{t) = Trse 



-f3H-{t) 



(12) 



where Tts denotes the trace over the system Hilbert 
space. Using Eqs. Q and pT|) . it follows that 



-fiH'it) 



r-i(i)TrBe-'5^(*) 



(13) 



where the right hand side coincides with the reduced den- 
sity matrix of the open system in thermal equilibrium 
with the heat bath. Again, t merely characterizes those 
parameter values that occur according to the protocol 
at time t. It does not indicate any dynamics. The ac- 
tual time dependent density matrix at time t does not, 
in general, coincide with e~^^ ^^'^ /Zs{t). 

We note that the Hamiltonian of mean force H* {t) typ- 
ically is a complicated operator-valued function not only 
of the system's parameters but also of the system-bath 
coupling strength, the bath temperature and possibly of 
other bath parameters. In the case of weak coupling the 
contributions from the bath and the interaction(s) are 
negligible and the Hamiltonian of mean force reduces to 
the bare system Hamiltonian [Oj] . 

In the classical limit, the quantum Hamiltonian of 
mean force becomes the classical Hamiltonian of mean 
force in Eq. This classical expression can further 

be simplified for a bath Hamiltonian consisting of a sum 
of potential and kinetic energies where the latter do not 
depend on positions and for an interaction that is inde- 
pendent of the bath momenta. Then the mere integration 
over the momenta yields identical factors in the numer- 
ator and denominator of Eq. pi[) . which cancel each 
other. The remaining term under the logarithm is then 
only a function of the system positions. This leads to the 
rcnormalization of the potential, - i.e. to the potential 



V*{x-t) 



Vsix,t) 
1 
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In 



J dyeM-PiVsiy) + VsB{x,y))) 



Jdyexp{-pVB{y)) 



(14) 



whereas the kinetic energy of the system remains un- 
changed by this procedure. We note that the given con- 
ditions are sufficient but not necessary in order that the 
potential of mean force captures the complete effect of a 
bath on the equilibrium properties of the system. Further 
simplifications result for example for a bath consisting of 
a set of harmonic bath oscillators which linearly couple 
to phase space functions of the system. Then, classically, 
the Hamiltonian of mean force coincides with the bare 
Hamiltonian of the system. This is in strong contrast to 
the behavior of quantum systems which couple as above, 
i.e.; linearly, to a harmonic heat bath. In this case, the 
Hamiltonian operator of mean force deviates from the 
bare system Hamiltonian with respect to the kinetic and 
the potential energy 

Conclusions.— Surprisingly enough, the Crooks theo- 
rem and the Jarzynski equality are valid for open quan- 
tum systems irrespectively of the coupling strength to 
their thermal environment and of the particular nature 
of their environment. These theorems, hence, are valid 
for all types of processes in which a classical or quantum 
system in contact with a thermal heat bath is driven 
out of equilibrium by classical, generally time-dependent 
forces. 
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